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Abstract 

This paper demonstrates the need for a nonconstant initialization for the fractional 
calculus and establishes a basic definition set for the initialized fractional differintegral. This 
definition set allows the formalization of an initialized fractional calculus. Two basis calculi are 
considered; the Riemann-Liouville and the Griinwald fractional calculi. Two forms of 
initialization, terminal and side are developed. 


1. Introduction 

The issue of initialization has been an essentially unrecognized problem in the 
development of the fractional calculus. Liouville’s definition ([5], p.21) for the fractional integral 
with lower limit of -°°and Reimann’s ([5], p.21) choice of the lower limit of c, were in fact 
related to the issue of initialization. Ross [8], [9] recognizes that to satisfy composition of the 
fractional differintegral, that the integrated function and its (integer order) derivatives must be 
zero for times up to and including the start of fractional differintegration. Podlubny ([7] pp. 
125,133) recognizes the need for initialization but carries it no further. Ross provides a history of 
the fractional calculus (see [6], p8.) in which he quotes A. Cayley referring to Riemann’s paper 
“The greatest difficulty in Riemann’s theory, it appears to me, is the interpretation of the 

complementary function ” Ross continues, “The question of the existence of a complementary 

function caused much confusion. Liouville and Peacock were led into error, and Riemann 
became inextricably entangled in his concept of a complementary function.” In retrospect the 
difficulties of Riemann over the role of the complementary function, which has been abandoned 
in this mathematics, may in fact have been related to the issue of initialization. The 
complementary function issue is raised here because an initialization function, which accounts for 
the effect of history, for fractional integrals and derivatives, will appear in the definitions 
presented. Its form is similar to Riemann’s complementary function, however, the meaning and 
use of this function is different (now clear). 

In the solution of fractional differential equations with an assumed history, it has been 
implicitly inferred by many authors ([1], [5], [6], [7], [10], and others), that an initializing 
constant, or set of constants, representing the value(s) of the fractional differintegrals (at t = 0) 
will provide an adequate representation for the effects of the past for each differintegral. That this 
is not true will be demonstrated in this paper. 

Finally, because the constant initialization of the past is insufficiently general, the widely 
used contemporary equation for the Laplace transform for the differintegral ([6], p. 135 for 
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example), based on that assumption also lacks sufficient generality. An alternative generalized 
form is presented. 

1.1 Proof of Non-Constant Initialization 

Consider the following qlh order fractional integrals of f (t ) . the first starting at time t = a , and 
the second starting at time t = c> a 

(')= (' - *r / W dx (11.1) 

and 

c (i f - r)‘ rl f(x) dx. (l . 1 .2) 

Assume that f(t) is zero for all t < a , then the time period between t = a and 
t = c maybe considered to be the “history” of the fractional integral starting at t = c , namely, 
c Then, we should expect that when this integral, { c d ~ q /(/)), properly initialized that 
it should function as a continuation of the integral starting at t — a . To achieve this an 
initialization must be addended to c d~ q f(t) so that the resulting fractional integration starting at 
t = c should be identical to the result starting at t = a for t > c . Thus, calling \jf the unknown 
initialization we have that 


c d; q f{t)+\if= a d; q f{t) t>c. 

(1.1.3) 

w= a d; q f{t)- c d; q f(t) t>c. 

(1.1.4) 

v= r ^)]( ? t)*7(t )dx- a d c q f{t) t> c . 

(1.1.5) 


Here y/ is seen to be a function of the independent variable t , completing the proof. We see that 
Iff is a generalization of the case for the ordinary integral ( q = 1 ), where 

c 

y/= jf (r)dr = constant. 


The insight for this behavior was originally obtained through the study of one- 
dimensional semi-infinite diffusion and wave equations [2], 

Having now recognized the need for a more general initialization, it must be decided if it 
is prudent to proceed as is done in the ordinary (integer order) calculus. That is, to append the 
initialization (constant of integration or constant initialization terms of ordinary differential 
equations) when required in an ad hoc manner or to formalize the process. 

Because of the increased complexity of the initialization relative to the integer order 
calculus case it is prudent to formalize the initialization, that is to include an initialization term 
into the definitions for the fundamental operators. The remainder of this paper will establish the 
definition set for an initialized fractional calculus, consider briefly the ‘Ross criteria” [8], [9] for a 
calculus, present a generalized (corrected) form for the Laplace transform of a differintegral, and 
demonstrate the solution of properly initialized fractional differential equations. 
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2. Definitions for the Initialized Fractional Calculus 

Several bases are possible for the initialized fractional calculus, these include the Riemann- 
Liouville, and the Griinwald formulations. The Riemann-Liouville form of the fractional calculus 
will be the only basis considered here. A consideration of a Griinwald based initialized fractional 
calculus may be found in [3]. In the development that follows, attention is restricted to real 
values of the order, q , of the various differintegrals. 

Two types of initializations are considered, ‘terminal initialization”, where it is assumed that 
the differintegral operator can only be initialized ("charged") by effectively differintegrating prior 
to the "start" time, t = c , and “side initialization”, where a fully arbitrary initialization may be 
applied to the differintegral operator at time t = c . The terminologies “terminal charging” and 
“side charging” have also been used to describe these initialization processes [3]. For discussion 
purposes, it is assumed that t (time) is the independent variable associated with the fractional 
differintegration and the function to be differintegrated is /(f). 

2.1 Initialization of Fractional Integrals 

2 . 1 . 1 T erminal I nitialization 

Terminal initialization is considered first. It is assumed that the fractional integration of interest 
"starts" at t = c (i.e. point of initialization). Further, /(f) = 0 for all t <a , and the fractional 
integration takes place for t>c>a. The initialization period (or space) is defined as the region 

a<t<c. 

The standard (contemporary) definition of a fractional integral will be accepted only when 
the differintegrand /(f) = Ofor all t < a 
Then, 

a D7 f('-rr/(r)</T, v>0 ,t>a, (2.1. Ll) 

r v / fl 

subject to /(f) = 0 for all f< a. The following definition of fractional integration will apply 
generally (i.e., at any f > c ) : 

c>aan d/(r)=0 \/t<a. 

The function y/(f,-v, a, c, t ) is called the initialization function and will be chosen such that 

aD?'f{t)= c D 7 /{*) t>c- (2. 1.1. 3) 

Substituting the results from equations 2. 1.1. land 2. 1.1. 2, then, for t > c and v ^ 0 gives 

-4-t f (t - r)'" l f(r)d t = -pr j(r - r)'" 1 /( r) d t+y/{f-v,a,c,t) (2.1.1 .4) 

r {v)i r (v)J 

r c r 

Because Jg(r)t/r= Jg(r)t/r + j g(z)di , 

a a c 

y/(f-v,a,c,t) = a D~' f(t) = f{z)dz , t>c, v>0. (2.1. 1.5) 

r ( v )« 


) ' f(z]dz + if/{f,-v,a,c,t ) t v>0,f > c, (2. 1.1. 2) 
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This expression for y/(t) applies for the terminal initialization condition. Clearly, If/ brings to 
the definition of the fractional integral the effect of the past, namely the effect of fractionally 
integrating /(f) from a to c . This effect will, of course, influence behavior after the time, 
t = c . The if/ function has the effect of allowing the function /(f) and its derivatives to start at 
a value other than zero, namely the value „ D~ q /(f) |,. f , and the I// -function continues to 

contribute to the differintegral response after f = c . That is, a function of time is added to the 
uninitialized integral, not just a constant. 

The integer order integrals under terminal initialization are of special interest. Evaluating 
equation 2. 1.1.5, for example, for v = 1 , indicates that !//(/ 1, o, c, f ) = constant. The general 
case is readily shown to be 

n 

\ f/(f,-n,a,c,t) = ^ l c i t' , n = 1,2,3- •• . 

/=0 

This, of course, is the same effect as seen in the integer order calculus using the '‘constant of 
integration.” It is important to note that the initialization of the qth fractional integral of f(t) is 
not unique in the following sense. That is, fit) can be considered as a composite function, for 
example, f(t)= g(t)U(t -c) + h(t)(U(t -ci)-U(t-c)), where £/(f)is the unit step function 

JO t < 0 

[l r>0. 

Then for this composite function /(f), it is the function g(t)U(t — c) that is being 
differintegrated and h(t)(U (t — a) — U{t — c)) is the function on which the initialization is based. 
This is analogous to choosing an arbitrary constant value to initialize (the integration of) dy / dt 
in the solution of an ordinary differential equation. 

2.1.2 Side Initialization 

When side initialization is in effect as opposed to terminal initialization equation 2.1. 1.2 is still 
taken as the operative definition and, 

V = v(0( i .e., is arbitrary) (2. 1 .2, l) 

That is, equation 2. 1.1.5 no longer (generally) applies. 



2.2 Initialization of Fractional Derivatives 

To extend the definition to the fractional derivatives, some issues must be addressed. The 
definition of the fractional derivative raises the following important questions in the context of 
initialization. Do fractional derivatives require an initializing function in general ? Further ; do 
integer order derivatives in this context require initialization functions? Clearly as we commonly 
think of derivatives, in the integer order calculus, the derivative is a local property and is 
represented geometrically as the slope of the function being considered and as such it requires no 
initialization. In the solution of differential equations the initialization constants which set the 
initial values of the derivatives really have the effect of accounting for the integration of the 
derivative from minus infinity to the starting time of the integration (of the differential equation). 

A study of the representation of semi-infinite systems using fractional differential equations 
[2], [3] indicates that the fractional derivative is not a local property as appears to be the case for 
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integer order derivatives (in the integer order calculus). Further, to solve such fractional 
differential equations an initialization function is required to handle the effect of the distributed 
initialization. Also, the integer order derivative in the fractional context can be shown to require 
an initialization function. Therefore, the answer to both questions is clearly YES. The impact of 
this is to require an initialization function for the fractional (and integer order) derivatives. 

Thus, a generalized integer order differentiation is defined as 

iffl , 

t>c, (2.2.l) 

dt 

where in is a positive integer and where !//(/, wi, a, c, f)is an initialization function. This is, of 

course, a generalization of the definition of the derivative, and for many cases, for example, 
usually in the integer order calculus, will be taken to be zero. It will be shown later that, for 

m = 1 with the condition of terminal initialization that t//(/, l,a,c,/) = 0. 

Now the uninitialized fractional derivative is defined as 

a D? f(t)^ a D™ a D; p f(f) <7 ^0,r>a,and /(f)=0Vf <a, (2.2.2) 

and (for convenience) m is the least integer greater than q, and q = m- p. Now as in the 
fractional integral case y/(f -p,a,a,t) = 0. Further, since lff(h,m,a,a,t) = 0, where 
h(t)= a D~ p f(t) , this definition specializes to the contemporary definition of the fractional 
derivative. 

Now the initialized fractional derivative is defined as 

c Dlfit)= c D7 f(t) cD? fit) 0, t>c>a, (2.2.3) 

where (for convenience) m is the least positive integer greater than q and q-m- p . 

The above definitions, equations (2.2.1) to (2.2.3), hold for both terminal initialization as well as 
side initialization. 


2.2.1 Terminal Initialization 

The initialization of the fractional derivative proceeds in a manner similar to the fractional 
integral, that is, it shall be required that; 

c D?f(t) = a D?f(t) Vt>c>a with ^>0. (2.2.1. l) 

Specifically, this requires compatibility of the derivatives starting at t — a and t = c, for t > c. 

It follows then that 

c D? c D; p f{t) = a D? a D; p f{t), q>0, t>c>a. (2.2.1.2) 

Expanding the generalized integral terms 



tY ' 1 f(T)dz + y{f-p,ci,c,t) 


( i f 

a D" -r--Ut-T) p ~ l f(T)dT+\l/{f-p,a,a,t] 

np)t 


(2.2. 1.3) 


Since ip ( f-p , a, aj) = 0 and using the definition (eq.(2.2. 1)) for the integer order derivative, 
equation (2.2. 1.3) may be written as 
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d" 

dt" 


1 ' 1 

' f( T )d T+ v(f-P’ a ’C, T )\ + ys(hi,m,a,c,t) = 

| f 

{ > c , (2.2.1.4) 


£_ 

dt m 


where h 1 = a D l p f(t) and h 2 = n D t p fit). The integer derivative is uninitialized at t = a, 

therefore yf[h z ,m,a,a,t) = 0. Then combining integrals gives 

( « c A 


d m ( 1 c 

vKfy , m, a, c, t) = — J (/ - T ) p ' [ f{x)dT - y/{f-p, a,c,t ) 


t > c. 


(2.2. 1.5) 


Under the condition of terminal initialization of the fractional integral, the argument of the 
derivative above is zero thus \ffi K h x ,m,a,c,t) = 0. For the case of side initialization of the 

fractional integral part of the fractional derivative y/(f,-p,a,c,t ) is arbitrary. Thus it can be 

seen from the above equation that either \f/(f ,-p,a,c,t) or y/{h l ,m,a,c,t) can be arbitrary but 
not both while still satisfying the requirements of the initialization (equation (2.2. 1.5)). 


2.2.2 Side Initialization 

The fractional derivative, side initialization case, can now be stated as 

,D,V(<)= r O;|^j(»-r) fH /(rVr+y(/,- ; >,a,c,r)|, q>0,t>c, (2.2.2.1) 

and m is the least positive integer >q with q = m- p, or equivalently as 

j m I t jm 

cD?f(t)= — /( r V T+ ^7^(/ -P,a,c,tyv(h,m,a,c,t\q>0,(2.2.2.2) 

where m is as above, t > c,and h(t) = a D~ p f(t) . Here both initialization terms are arbitrary 
and thus may be considered as a single (arbitrary) term, namely 

d m 

y/{f,q,a,c,t)=—-y/(f,-p,a,c,t)+y/(h,m,a,c,t). ( 2.2.23 ) 

d t 


In the case of terminal initialization of the fractional integral part of equation (2.2.2. 1) 

\j/{f ,-p,a,cj) will be as defined by equation (2. 1.1.5). It is noted, that the a and c arguments 
in the y/ functions in equations (2.2.2. 1) to (2,2.2. 3) are carried for parallelism to the previous 
cases and are not intended to infer an initialization period under side initialization. 


2.3 A Simple Example 

A simple example will be helpful. Consider the semi-integral of /(f) = (f — a)U(t-a), 
then from ([6], pp. 63-64) the uninitialized semi-integral becomes 

a D; 1/2 (t-ayj{t-a)= , t>a. (2.3.l) 

Now, initializing from the point t — c 
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Semi-Integrals & Initialization Function 


t>c 


= l>c. (2.3.2) 

Consider now the terminal initialization, 

V<(/-l/2,a,c,i)=-p^y^j(r-r)"'' : (r-aVr i>c, 

= ^j)[(>-cr(- 2 ' + 3«-^)-('-«) l ' ! (- 2 ' + 2‘')]. '>* (2-3.3) 

The numerical evaluations of these equations for specific numerical values, a = -1, c = 1, are 
shown in the graphs of figure 1 . 



-1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 3 3.5 


Time 

Figure 1 Sample Problem — Semi-Integral of (t-a) vs Time with 

a= -1, c=l 

3. Criteria for A Fractional Calculus 

Ross [9] provides a set of criteria for a fractional calculus. The criteria include the 
following properties: backward compatibility with the normal (integer order) calculus, the zero 

property, namely c D l °/(f)= f(t\ linearity of the operators, and that the index law should hold, 

that is, C D“ c DJf(t)= c DJ c D‘‘f(t)= c D' + ''f(t). Under the conditions of terminal initialization 
the above properties are each shown to hold ([3], pp. 26-44). This provides credibility to the 
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initialized fractional calculus, also completion of the proofs yields constraints on the allowable 
initialization functions. For example, linearity of the fractional integral 

c D; r (bf(t)+k g (t))=b c D; v f(t)+k c D;' g (t) t>c, ( 3 . 1 ) 

holds if and only if 

yf {bf +kg-\\a,c,t)= by/(f -v,a,c,t)+kyr(g-v,a,c,t) t>c. (3.2) 

4. Laplace Transform of Fractional Differintegrals 

The Laplace transform of the differintegral is an important tool in the solution of 
fractional differential equations. The following form for the Laplace transform of the fractional 
derivative is given by many authors (see for example, ([7], p.105), ([6], p. 134), or an equivalent 
form in ([5], p. 1 23)) 

n-l< p< n . (4.1) 

k= 0 

The Laplace transform of the fractional integral is given ([7], p.104), and ([6], p.134) as, 
L{ D D; r ‘f(t)}=s- p F(s) p> 0. (4.2) 

The form of the summation in equation (4.1) infers that the p- th order derivative is being 
decomposed into a fractional integral (order < 1) and (n-1) order 1 derivatives. Further inferred is 
that each derivative is initialized by an impulse at t = 0 . This situation is a residual from the 
integer order calculus and lacks sufficient generality for a properly initialized fractional calculus. 

It should be noted that in equations (4.1) and (4.2) above D p and D~ p refer to the contemporary 
uninitialized fractional derivative and integral respectively. 

For the initialized fractional calculus the Laplace transform of the initialized fractional 
differintegral ([3], p.61) is given by 

LlD*f{t)}=s«L{f{t)}+Lfy(j,q,aAt)\ Vq, (4.3) 

where y/ depends on; the assumed past history of 0 D? f(t), and the assumed decomposition of 

0 D‘1 f (t ) (as required in the problem definition). It should be noted that for q > 0 in the most 
general case (side initialization) that the equivalent form for y/ must be used. That is 

im 

y/(f,q,a,0j)=-—yf(f-p,a,0,t)+y/(h,m,a,0,t) q> 0, where h(t)= a d; p f(t). (4.4) 
a t 

More powerful forms than equation (4.3) have been derived that account for all possible 
decompositions of 0 D, f (f) . For a detailed explanation of the decomposition issue the reader is 
referred to ([3], p.46-58). 


5. Fractional Differential Equations 

Proper initialization is crucial in the solution and understanding of fractional differential 
equations. The application of the initialized fractional calculus to the solution of initialized 
fractional differential equations will be illustrated with the following examples. 
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5.1 Example 1 

Podlubny ([7], p. 138) and Oldham and Spanier ([6], p. 157) consider the following 
fractional differential equation 

0 d;'7W+*>/(')=o, / > o; O o; , '7(d«=c. (s.i.i) 

The notation 0 D)' 2 fit ) is that of Podlubny, and of course, refers to the uninitialized derivative. 
Equation (4. 1) above, is used by Podluny, to obtain the Laplace transform of equation (5. 1 . 1) as 


Ftt-TBT 


(5.1.2) 


s'" +b 

The inverse transform is then given in terms of a two-parameter Mittag-Leffler series expansion 


[L (-byfi). 

>’ 2 

(5.1.3) 


(5.1.4) 


For b = 1 , then this is determined to be 


which agrees with the result of Oldham and Spanier [6]. This is contrasted with the following 
approach using the results from the initialized fractional calculus. 

We now solve equation (5.1.1) again, but now 0 D)' 2 fit ) is interpreted as an initialized 
fractional derivative as defined in part 2.2 above. Thus, we have 

0 Df 2 f{t)+bf(t)=0, t >0, i/s(f,l/2,a,0,t) is arbitrary. (5.1.5) 

This may be rewritten as 

0 d l / 2 f(t)+y/(f,l/2,a,0,t)+bf{t)=0, t> 0, y/(/,l/2,a,0,f)is arbitrary. (5.1.6) 

The Laplace transform of equation (5.1.5) using equation (4.3) is 

s U2 +b s U2 +b' } 

This equation should be contrasted to equation (5.1.2) above, they are only the same when 
y/[t) = -C<5(V), that is, when an impulse at t = 0 is used to initialize the fractional differential 
equation! 

Now the /^-function and its Laplace transform [4], are given by 

» / \n f \{n+\)q-v-\ v 

7? ,(a,c,t)=^^y — — r- o -f — , Re(^-v)>0, Re(5)>0. (5.1.8) 


The general inverse for equation (5.1.7), is obtained by applying Laplace convolution integral, 

t 

fit) = -jR U20 {-b,0,t-T)y/(r)dT t> 0. (5.1.9) 

0 

Thus, with arbitrary l/r(f)this provides the most general solution to equation (5.1.1), or (5. 1.5). If 
we take \(/{t)= —Cd(t) inequation (5.1.9), the result is 

f(t)=CR u20 {-b, 0,/), (5.1.10) 

which is identical with the result of equation (5.1.2). In the context of the initialized fractional 
calculus, this might be considered as a pathological result of little interest. A more useful result 
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would be obtained under the assumption of terminal initialization, that is apply equation (4.4) 
with q-m-p , m = 1, p = 1/2, and \j/(h,m,a,0,t) = 0 . Then 

V/(/,l/2, a ,0,r)=|;y(/,-l/2,a,0,r)=|.y^j(r-rr ! /(rVr. (51. ll) 

In this evaluation ( /(f), t < 0) need not be identical to ( f(t ), t > 0), if f{t ) is considered to be a 
composite function. 

5.2 Example 2 

A second fractional differential equation is considered by Podlubny ([7], p. 139), and has 
a special case considered by Oldham and Spanier ([6], p. 159 Q = 1, q = 1 / 2 ). The general case, 
of which it is said “encounters great difficulties except when q — Q is integer or half-integer,” is 
given by 

0 Df/(/)f 0 O,V(')=*('). c = [ 0 D,°-7('KDr7WL. (5.2.1) 

Again, the notation 0 D, 2 fit) is that of Podlubny, and refers to the uninitialized derivative. Here, 
we shall assume the notation is that of this paper, that is represents the initialized fractional 
derivative, and we shall relax the requirement that q-Q '\s integer or half-integer. The 
initialization of equation (5.2.1) will be replaced by two separate ones to identify the most general 
solution. Then the Laplace transform of equation (5.2.1) gives 

s Q F(s)+s 9 F(s) = H{s)-y/ l (f,Q,a,0,s)-\i/ 2 (f,q,a,0,s), (5.2.2) 

where the t//'s are subscripted for convenience. This may be written, assuming that Q> q, as 

F(s)= -J— (5.2.3) 

The solution is easily written using the /^-function as, 

= ] -TXhiryy/Xf ,Q,a,0,T)-y/ 2 (f ,q,a,0,r))dT . (5.2.4) 

0 

This solution may be specialized to that of Podlubny by taking y/ l = — C t S(t) , \ j/ 2 = —C z 5(t ) 
and C = C, + C, , namely 

t 

/(l)=-a H .,(-l,0,l)+ 1 (-1,0, l-z)(h(t))dr . (5.2.5) 

0 

The important issue here is that the form of equation (5.2.5) does not allow the effect of 
continuing the past as does equation (5.2.4). 

5.3 Example 3 

In this example, an approach will be demonstrated in which the entire fractional 
differential equation is initialized as opposed to the above examples where the individual 
derivatives were initialized. We now generalize equation (5.2.1) as follows 

(5.3.1) 

This may be rewritten as 

c d?f(t)+ c d?f(t)= h(t)-y/Xf,Q,a,c,t)-y z (f,q,a,cs)=h{t)-y eqiv {^ . (5.3.2) 

Now the Laplace transform of the uninitialized fractional derivative is ([3], p.60) 

lU 5 /(/)}= lU'("M)/('))}= f-vq/(r+fc)} . (5.3.3) 


NASA/TP— 2000-209943 


10 



Thus, the Laplace transform of equation (5.3.2), for Q> q, is given by 


M/(f+c)}- = ‘■"G(sih(s)-V, t ,(s)) (5.3.4) 

^ s Is I J. I 

Clearly the behavior of the fractional differential equation (or the system it represents) is captured 
by the function G(s) = l/(s 9 (s' 12 * 9 + 1)). We now consider two time domain segments, for domain 
(1) we take a, = c, = 0 , and for domain (2) a 2 = 0, c 2 = 1 . 

Domain (1) 

Now domain (1) will be used as the initializing period for domain (2). To do this let 
— l)), for simplicity here we take h*(t)=l . Wealsotake 
/(f) = 0, Vf <0thereby inf erring t/r ] (/,g,0,0,r)= y/ 2 (f ,q,0,0,t)= 0. Thedomain(l) 
solution then is given by 


L{fM- 


1 

1 

1 

1 

1 

?+ V-«+i) 

l * ) 


to. 


(5.3.5) 


and 

/, (0 ! (- 1A /)- u(f - 1 )R Q . q .. q . l (- 1 A t - 1) t> 0. (5.3.6) 

For 0 < t < 1, then /, (/) = R Q . q . q . x (- 1,0, t). 


Domain (2) 

Now for domain (2), t > 1, we consider the initialization period to be, domain (1) above, 
0 < t < 1, therefore, a 2 = 0, c 2 = 1 , and again we take f(t) = 0 , V t < 0 . Further, for clarity, 
we only consider the unforced problem, that is we take h 2 (t )= 0, therefore from equation (5.3.2) 
we have 

^f 2 (t)+^f 2 {t) = -^(/ 2 ,GAl,0-^(/2^Al,0 = -^(0- (5-3.7) 

The Laplace transform, using equation (5.3.3), is given by 

-vyCO 
Tv^+I) 

The initialization If/ v (t) is now chosen based on /t, (0 < t < 1) , that is the historic forcing 


+ !)} = • 


(5.3.8) 


function, thus 

Weqiv = = l)), 

hence 



Substituting this result into equation (5.3.8) 



(5.3.9) 

(5.3.10) 


*-■ /-{/=(<+ 1)}= Jr I'a = — gw = a{/;M}. (5.3.1 0 

S 15 t 1 1 S 

Applying the Laplace shifting theorem, we have 

(5-3.12) 

Thus for t > 1 we have the important result /,(?)= /, (f ) t > 1 , as would be expected from a 
proper initialization theory. It is not possible to obtain such results from the contemporary theory 
referenced earlier. 
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Summary 


This paper proves the need for an initialized fractional calculus. The paper presents the 
definition sets required for initialized fractional calculi. Two underlying bases have been used, 
the Riemann-Liouville based fractional calculus and the Grunwald based fractional calculus (by 
reference). 

The significant result is that when fractionally differintegrating with respect to t that a 
function of t is required as an initialization as opposed to the constant initialization used in the 
integer order calculus. 

Two types of initialization are introduced “terminal initialization” and “side 
initialization”. Proofs of Ross’ criteria and the initialization constraints that ensue from the 
criteria have been referenced. Corrected forms for the Laplace transforms for fractional 
differintegals that properly account for the initialization function have been presented. 
Commonly studied fractional differential equations have been solved to demonstrate the various 
aspects of initialization. 
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